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Abstract 

Background: The evolutionary history of the biota of North Africa and Arabia is inextricably tied to the complex 
geological and climatic evolution that gave rise to the prevalent deserts of these areas. Reptiles constitute an 
exemplary group in the study of the arid environments with numerous well-adapted members, while recent studies 
using reptiles as models have unveiled interesting biogeographical and diversification patterns. In this study, we 
include 207 specimens belonging to all 12 recognized species of the genus Stenodactylus. Molecular phylogenies 
inferred using two mitochondrial (12S rRNA and 16S rRNA) and two nuclear {c-mos and RAG-2) markers are 
employed to obtain a robust time-calibrated phylogeny, as the base to investigate the inter- and intraspecific 
relationships and to elucidate the biogeographical history of Stenodactylus, a genus with a large distribution range 
including the arid and hyper-arid areas of North Africa and Arabia. 

Results: The phylogenetic analyses of molecular data reveal the existence of three major clades within the genus 
Stenodactylus, which is supported by previous studies based on morphology. Estimated divergence times between 
clades and sub-clades are shown to correlate with major geological events of the region, the most important of 
which is the opening of the Red Sea, while climatic instability in the Miocene is hypothesized to have triggered 
diversification. High genetic variability is observed in some species, suggesting the existence of some undescribed 
species. The 5. petrii - S. stenurus species complex is in need of a thorough taxonomic revision. New data is 
presented on the distribution of the sister species S. sthenodactylus and S. mauritanicus. 

Conclusions: The phylogenetic hypothesis for the genus Stenodactylus presented in this work permits the 
reconstruction of the biogeographical history of these common desert dwellers and confirms the importance of the 
opening of the Red Sea and the climatic oscillations of the Miocene as major factors in the diversification of the 
biota of North Africa and Arabia. Moreover, this study traces the evolution of this widely distributed and highly 
specialized group, investigates the patterns of its high intraspecific diversity and elucidates its systematics. 
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Background 

North Africa and Arabia are home to a unique fauna 
and flora that has been shaped by the combination of 
several factors including the harsh climatic conditions 
of the Sahara and Arabian deserts, the episodic appea- 
rance of humid cycles, and by the complex geological 
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evolution of the area [1-9]. One of the most important 
geological phenomena of the entire Cenozoic that oc- 
curred in this area was the break-up of the Arabian plate 
from Africa. Tectonic activity started approximately 
30 Ma ago at the central Gulf of Aden with the forma- 
tion of a rift basin in the Eritrean Red Sea and initial rif- 
ting at the Afar zone. A second phase of volcanism 
occurred 24 Ma ago, causing extension and rifting 
throughout the entire Red Sea, from Yemen to Egypt, 
as well as uplifting of the newly-formed continental 
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shoulders [1]. Nevertheless, fluctuations of the sea level 
during the Miocene permitted the formation of transient 
land connections [1,10] that were subsequently lost [11]. 

The establishment of the Afro-Arabia - Eurasia land 
bridge {Gomphotherium bridge) was another crucial 
event with major biogeographical implications [12-14]. 
Following the opening of the Gulf of Aden and the Red 
Sea and with the counterclockwise rotation of the Arabian 
plate, a first connection was presumably formed between 
the latter and the Anatolian plate, and subsequently 
with Eurasia. Although the connection between the 
Mediterranean Sea and the Indian Ocean is hypothe- 
sized to have been re-established in the Upper Middle 
Miocene, around 15 Ma ago, it is believed that posterior 
to this date the land bridge has been continuously present 
[15]. Important faunal and floral exchanges have been 
attributed to the establishment of this connection ([12-14] 
and references therein). 

Although the origin of the Sahara and Arabian deserts 
is still hotly debated [16-19], it is generally accepted that 
climatic development in the late Miocene, as a result of 
major growth of the East Antarctic Ice Sheet and polar 
cooling, lead to an increase in aridification of mid- 
latitude continental regions [4] and that this had a pro- 
found effect on the diversification of faunas [20-22]. 

Reptiles are among the commonest inhabitants of arid 
areas and have long been used in biogeographic, ecological 
and evolutionary studies [23], constituting thus excellent 
models to investigate how diversity is originated and main- 
tained. Several cases of faunal exchanges in both directions 
between North Africa and Arabia have been described 
(e.g. [2,13,24]) showing that there is not a single pattern, 
but rather different hypotheses including both vicariance 
and dispersal, heavily dependent on the estimated time- 
frame of the events. Moreover, several studies have shown 
that climatic changes towards aridity and contraction/ 
expansion of the Sahara and Arabian deserts have played 
a decisive role in reptile species diversification [25-29]. 

Gekkonid lizards of the genus Stenodactylus Fitzinger, 
1826 [30] are one of the most characteristic and abun- 
dant elements of the fauna of the arid and hyper-arid 
regions of Arabia and North Africa [31]. The genus 
comprises twelve species that are distributed in a more 
or less continuous range across northern Africa and 
Arabia, with an apparently isolated population in north- 
ern Kenya and extending around the Arabian Gulf to 
coastal southwestern Iran ([32,33]; see Figure 1). Up to 
three species may occur at a single locality and, where 
such sympatry exists, resource partitioning is largely 
achieved by microhabitat segregation, with species occu- 
pying different soil types [34]. Gravel plains, hard sand 
and aeolian soft sand all have their characteristic species 
that show specialized morphological adaptations. These 
include the presence of depressed and fringed toes, which 



increase the surface area and improve grip in the aeolian 
sand dune specialists Stenodactylus doriae (Blanford, 1874 
[35]), S. petrii Anderson, 1896 [36] and S. arabicus (Haas, 
1957 [37]). Extensive webbing is also observed between 
the fingers for efficient sand burrowing in S. arabicus 
[31,32,38]. When two species are regularly found on the 
same substrate, they greatly differ in size and there are 
corresponding differences in the size of prey taken [32] . 

Morphologically, Stenodactylus is fairly homogeneous 
and all species exhibit phalangeal reduction that pro- 
duces a formula of 2.3.3.4.3 on both fore and hind limbs 
and are also characterized by a very high scleral ossicle 
number (20-28) [31,39]. A morphology-based phylogen- 
etic hypothesis has been proposed by Arnold (1980) 
[31]. Although these two characters are also present in 
Pseudocemmodactylus khobarensis Haas, 1957 [37], which 
was widely accepted as a Stenodactylus member [31,39], a 
recent phylogenetic study by Fujita and Papenfuss (2011) 
[40] including specimens of the former and six out of the 
twelve species of the genus Stenodactylus proposed the 
resurrection of the genus Pseudoceramodactylus. This was 
done in order to deal with the resulting paraphyly of 
Stenodactylus, caused by the branching of two represen- 
tatives of the genus Tropiocolotes between P. khobarensis 
and the six Stenodactylus included in their analyses. Their 
molecular analyses also uncovered high levels of genetic di- 
vergence between the different Stenodactylus species. Ge- 
netic variability within some of the species, like S. arabicus 
and S. doriae, was also high and this could be linked to 
biogeographic discontinuities among some of the hyper- 
arid areas in Arabia. 

Although Stenodactylus includes a relatively low number 
of species compared to other gecko groups in these areas, 
such as Pristurus, Tarentola or Hemidactylus [26,41-46], 
its relatively high level of resource partitioning and 
habitat specialization has allowed the different species 
to successfully colonize almost all available habitats in 
the arid and hyper-arid regions of North Africa and 
Arabia. It constitutes, therefore, a very interesting, but still 
poorly studied, genus that makes an excellent model for the 
study of desert biodiversity and biogeography. The main 
objectives of the present work are: (1) to provide for the 
first time a complete phylogeny of the genus Stenodactylus 
and evaluate its concordance with previous molecular 
and morphology-based studies; (2) to investigate the bio- 
geographical and diversification patterns of Stenodactylus; 
and (3) to explore the interspecific relationships, the pat- 
terns of intraspecific diversity and the possible presence of 
unrecognized divergent lineages in Stenodactylus. 

Methods 

Taxon sampling, DNA extraction and sequencing 

A total of 207 individuals of Stenodactylus representing 
all twelve currently recognized species were included in 
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Figure 1 Sampling localities of the Stenodactylus specimens used in this study. Colors and locality numbers refer to specimens in Figures 2 
and 3 (see also Additional file 1: Table SI). The global distribution of the genus is seen on the upper right (data from Sindaco and 
Jeremcenko, 2008). 



the present study. Whenever possible, we tried to in- 
clude multiple populations for each species in order to 
assess intraspecific variation; sampling was especially 
intense in the three African species with very large dis- 
tribution ranges. In addition, three Pseudoceramodactylus 
khobarensis and eight individuals representing six species 
of the genus Tropiocolotes were included in an attempt 
to further test the relationship between Stenodactylus, 
Pseudoceramodactylus and Tropiocolotes. Four additional 
specimens from other closely related genera [47-49] were 
used as outgroups and sixteen specimens, from several 
genera, were added in order to estimate divergence times 
(see below). Additional file 1: Table SI lists all 238 samples 
used in the present work with their extraction codes, vou- 
cher references, localities and GenBank accession num- 
bers (KC190516-KC191151). 

Genomic DNA was extracted from ethanol-preserved tis- 
sue samples using the DNeasy Blood & Tissue Kit (Qiagen, 
Valencia, CA, USA). All 222 specimens included in the 
phylogenetic analyses were sequenced for two mitochon- 
drial gene fragments: 378-388 base pairs (bp) of 12S rRNA 
(12S) and 498-536 bp of the 16S rRNA (16S). A subset of 
106 specimens, including representatives from all inde- 
pendent lineages recovered by the analysis of these two 



fragments, was also sequenced for two nuclear markers: 
660 bp of the oocyte maturation factor MOS (c-mos), and 
410 bp of the recombination activating gene 2 {RAG-2). 
Primers used for the amplification and sequencing of the 
12S, 16S, c-mos and RAG-2 gene fragments as well as PCR 
conditions applied in the present work are listed in detail in 
Table 1. All amplified fragments were sequenced for both 
strands. Contigs were assembled in Geneious v.5.3 [50]. 

Phylogenetic analyses and hypothesis testing 

DNA sequences were aligned using the online version of 
MAFFT v.6 [51] with default parameters (gap opening = 
1.53, offset value = 0.0) for the mitochondrial genes and 
with modified parameters (offset value = 0.1) for the nu- 
clear genes, in which long gaps are not expected. Coding 
gene fragments {c-mos and RAG-2) were translated into 
amino acids and no stop codons were observed. Uncor- 
rected ^-distances were calculated in MEGA v.5 [52]. 

Phylogenetic analyses of the combined dataset were 
done employing maximum likelihood (ML) and Bayesian 
(BI) methods. Separate ML analyses were also performed 
on 12S, 16S, c-mos and RAG-2 to test for conflicting sig- 
nal among genes. Best-fitting nucleotide substitution mo- 
dels were selected for each partition under the Akaike 
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Table 1 Primers used in this study 



Gene fragment 


Primer name 


Or. 


Sequence (5'- 3') 


Reference 


PCR conditions 








12S 


1 2Sa 


F 


AAACTGGGATTAGATACCCCACTAT 


Kocher et al. (1989) 


94° (5 1 ) 


■ 940 (45 


51° (45"), 72° 


(80") : 


X 35; 72 (5') 




12Sb 


R 


GAGGGTGACGGGCGGTGTGT 


Kocher et al. (1 989) 














L1.STEN0 


F 


GGATTAGATACCCCACTATGC 


This study 


94° (5 1 ) 


■ 940 (45") i 


52° (45")', 72 c 


(90") X 35; 72° (5') 




H1.STEN01 


R 


TGACGGGCGGTGTGTACG 


This study 












16S 


16Sa 


F 


CGCCTGTTTATCAAAAACAT 


Palumbi (1996) 


940 


■ 940 (45 »^ 


51 (45"), 72 (J 


50") X 35; 72° (5') 




1 65b 


R 


CCGGTCTGAACTCAGATCACGT 


Palumbi (1996) 














16SaST 


F 


ATCAAAAACATCGCCTTTAGC 


This study 


940 


■ 940 (45 »^ 


57° (45"), 72° 


(70") : 


X 35; 72° (5 1 ) 




16SbST 


R 


CTGAACTCAGATCACGTAGGAC 


This study 












C-mos 


rUr 


r 

r 


1 1 1 Ljb 1 1 Lr\o 1 L 1 ALAAbbL 1 AL 


bamule et al. [zUUoj 


940 (5.) 


; 94° (30"), 


55° (45"), 72° 


(70") : 


X 35; 72° (10') 




FUR 


R 


AGGGAAGATCCAAAGTCTCCAAT 


Gamble et al. (2008) 














G73_STENO 


F 


GCTGTAAAGCAGGTGAAGAAATGC 


This study 


94° (5') 


; 94° (45"), 


56° (45"), 72° 


(80") : 


X 35; 72° (5') 




G74_STENO 


R 


GAACATCCAAAGTGTCCAATC1TGC 


This study 














G73.5_STENO 


F 


GCATTTGGACTTAAAACCTG 


This study 














G708 


R 


GCTACATCAGCTCTCCARGA 


Hugall et al. (2008) 












RAG-2 


RAG2-PY1-F 


F 


CCCTGAGTTTGGATGCTGTACTT 


Gamble et al. (2008) 


940 (5.) 


; 94° (45"), 


55° (45"), 72° 


(70") : 


X 35; 72° (5') 




RAG2-PY1-R 


R 


AACTGCCTRTTGTCCCCTGGTAT 


Gamble et al. (2008) 













List of primers used in the amplification and sequencing of gene fragments, with the corresponding source and PCR conditions. 
1 0rientation. 



information criterion [53] using jModelTest v.0.1.1 [54]. 
The GTR + I + G model was independently estimated for 
each of the 12S, 16S, RAG-2 partitions and the GTR + G 
model for the c-mos partition. Alignment gaps were trea- 
ted as missing data and the nuclear gene sequences were 
not phased. Hemidactylus frenatus was used for rooting 
the tree, based on published evidence [47,48]. 

A Bayesian analysis of the combined dataset was per- 
formed in MrBayes 3.1.2 [55,56] with best fitting models 
applied to each partition (gene) and all parameters 
unlinked across partitions. Analyses ran for 10 7 genera- 
tions, with sampling intervals of 1000 generations, produ- 
cing 10000 trees. Convergence and appropriate sampling 
were confirmed examining the standard deviation of the 
split frequencies between the two simultaneous runs and 
the Potential Scale Reduction Factor (PSRF) diagnostic. 
Burn-in was performed discarding the first 2500 trees of 
each run and a majority- rule consensus tree was gener- 
ated from the remaining trees. ML analyses were per- 
formed in RAxML v.7.0.3 [57]. A GTR + I + G model was 
used and parameters were estimated independently for 
each partition. Node support was assessed by bootstrap 
analysis [58] including 1000 replications. 

Haplotype networks were constructed for the two nu- 
clear markers c-mos and RAG-2. The software PHASE 
v.2.1.1 [59,60] was used to resolve the haplotypes where 
more than one heterozygote position was present. Input 
files were prepared using Seqphase [61]. In order to in- 
clude as much information as possible for the better 
resolution of the haplotypes, the alignment of all full- 
length sequences of each marker was used. Phase 



probabilities parameter was set at 0.7 and all other set- 
tings were set by default. The network of the resulting 
haplotypes was calculated with TCS v. 1.21 [62] applying 
default settings (probability of parsimony cutoff: 95%). 

Topological constraints to test alternative topologies 
were constructed by hand and compared to the uncon- 
strained (best) tree using the Approximately-Unbiased 
(AU) [63] and Shimodaira-Hasegawa (SH) [64] tests. 
Per-site log likelihoods were estimated in RAxML 7.0.3 
[57] and P- values were calculated using CONSEL [65]. 
Tests were also run in a Bayesian framework, where the 
relative support of competing hypotheses given the data 
was quantified using the Bayes factor (BF) [66]. Topolo- 
gies were constrained in analyses run in BEAST v.1.6.1 
[67], the marginal likelihood for each topology was esti- 
mated using the harmonic mean estimator and the Bayes 
factors were calculated by taking the ratios, as estimated 
in Tracer v.1.5 [68]. 

Estimation of divergence times 

A Bayesian approach was used to estimate divergence 
times as implemented in the software BEAST v.1.6.1. 
The dataset comprised sequences from all four partitions 
(the nuclear genes c-mos and RAG-2 unphased). An ar- 
bitrarily pruned phylogeny was used in order to include 
only one representative from each species or main lineage 
uncovered with the concatenated analysis (45 specimens 
in total; see Additional file 1: Table SI). This method 
excludes closely related terminal taxa because the Yule 
tree prior does not include a model of coalescence, which 
can complicate rate estimation for closely related 



Metallinou et al. BMC Evolutionary Biology 2012, 12:258 
http://www.biomedcentral.com/1471-2148/12/258 



Page 5 of 1 7 



sequences [69]. Additionally, several taxa belonging to 
other gecko genera were added for the calibration process 
(see below). 

Two individual runs were performed for 4 x 10 7 genera- 
tions with a sampling frequency of 4000 and the results 
were combined to infer the ultrametric tree after discard- 
ing 10% of the samples from each run. Models and prior 
specifications applied were as follows (otherwise by de- 
fault): GTR + I + G (12S, 16S), GTR + I (c-mos), HKY + I 
{RAG-2); Relaxed Uncorrelated Lognormal Clock (esti- 
mate); Yule process of speciation; random starting tree; 
alpha Uniform (0, 10); yule.birthRate (0, 1000). Parameter 
values both for clock and substitution models were 
unlinked across partitions. 

Unfortunately, no fossils belonging to Stenodactylus, 
Pseudoceramodactylus or Tropiocolotes are known, pre- 
cluding the direct estimation of the time of the clado- 
genetic events within our study group. Consequently, 
the estimation was based on well-known calibration 
points published in recent literature [70,71] related to 
members of the families Phyllodactylidae and Sphaero- 
dactylidae (see Additional file 1: Table SI). Three fossil 
and biogeographical calibration points were applied as 
"soft" priors, in order to account for uncertainty in the 
date of the corresponding nodes: (1) the minimum age 
for the divergence between Euleptes and its sister clade 
was set to 22.5 Ma ago using the approximate age of a 
fossil Euleptes [72,73] (Lognormal distribution: median 
22.5, 97.5% 36.55); (2) the split between Teratoscincus 
scincus - Teratoscincus roborowskii caused by the Tien 
Shan-Pamir uplift 10 Ma ago [74-76] (Lognormal distribu- 
tion: median 10.08, 97.5% 12.96); (3) the age of El Hierro 
island [77] at 1.12 Ma ago, assuming that divergence be- 
tween Tarentola boettgeri hierrensis and Tarentola boettgeri 
bischqffi began soon after its appearance [26,44] (Uniform 
distribution: lower 1, upper 1.12). In order to cross-check 
the results, the posterior mean rates of the mitochondrial 
gene fragments of our analysis were compared to the 
rates calculated for well-known and well-studied reptile 
groups from the Canary Islands (the geckos of the genus 
Tarentola, the lacertid lizards of the endemic genus 
Gallotia and the skinks of the genus Chalcides), for which 
robust calibrated phylogenies have been produced in se- 
veral independent analyses ([26,45,78-80], among others), 
and evolutionary rates for the 12S gene have been 
obtained using BEAST [44]. 

Ancestral area reconstruction 

MacClade v. 4.08 [81] was used to reconstruct the ancestral 
areas for the Stenodactylus species in a parsimony frame- 
work, using both delayed transformation (DELTRAN) and 
accelerated transformation (ACCTRAN). Additionally, in 
order to incorporate branch-length information, ML was 
used as implemented in the Mesquite software package 



[82]. Both Markov k-state 1-parameter and Asymmetrical 
Markov k-state 2-parameter models were applied and a 
likelihood ratio test was used to choose the best reconstruc- 
tion. Two states, Arabia and Africa, were identified in the 
extant species depending on the present distribution of the 
species [33] and were used with both methodologies. 

Results 

Phylogenetic analyses and topological tests 

Two datasets were used to infer the phylogenetic rela- 
tionships of the genus Stenodactylus: a mitochondrial 
one for building the preliminary phylogeny and analy- 
zing the divergence patters, and a multi-locus one for 
producing a more robust phylogeny (TreeBASE ID: 
13567). The first dataset consisted of an alignment of 
974 bp of mitochondrial DNA (415 bp of 12S and 
559 bp of 16S, of which 270 in both cases were variable 
positions) for 222 terminals including 207 Stenodactylus. 
The results of the ML and BI of this dataset were very 
similar and are summarized in Supplementary Figure 1 
(Additional file 2: Figure SI). In order to improve our 
phylogenetic hypothesis applying a multi-locus approach, 
a second dataset was assembled with a selection of 
106 terminals, including 91 Stenodactylus (see Additional 
file 1: Table SI) for which two extra nuclear genes were 
sequenced. The aligned dataset consisted of 2092 bp 
(419 bp of 12S, 560 bp of 16S, 703 bp of c-mos and 
410 bp of RAG-2, of which 262, 269, 109 and 99 
positions were variable, respectively). The result of 
the phylogenetic analyses of the concatenated align- 
ment of four genes is shown in Figure 2. Well-supported 
relationships in the independent gene trees were 
congruent among partitions, but at this level not all 
markers offered sufficient resolution to differentiate par- 
ticularly between S. sthenodactylus and S. mauritanicus 
(data not shown but see below). The networks con- 
structed for the phased haplotypes of the nuclear 
markers are presented in Figure 3. Not all ambiguities 
were resolved. 

Both ML and Bayesian analyses of the concatenated 
alignment of four gene fragments (Figure 2) gave almost 
identical results to the mtDNA tree from the Additional 
file 2: Figure SI There is low support over the relationships 
between the genera Stenodactylus, Pseudoceramodactylus 
and Tropiocolotes. According to the results, the North 
African T. algericus and T. tripolitanus branch first and 
P. khobarensis is sister to a poorly supported clade 
formed by two reciprocally monophyletic groups: one in- 
cluding T. scorteccii, T. steudneri, T. nubicus and the Middle 
Eastern T. nattered and the other one including all 12 spe- 
cies of the genus Stenodactylus. In order to further investi- 
gate these relationships, three topological tests were carried 
out: (1) Stenodactylus + Pseudoceramodactylus were forced 
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i S. pulcher - 1 - Yemen [47] 
i* " - 




Crossobamon orientalis 
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S. pulcher - 2 - Yemen [48] 
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S. arabicus - 6 i r <n 
-S. arabicus -5 I 0ma "I8] 
S. leptocosymbotes - 2 - Oman [21] 
S. leptocosymbotes - 3 - Oman [22] f 
S. leptocosymbotes - 4 - Oman [23] 1 
S. leptocosymbotes ■ 1 - Oman [20] 
S. leptocosymbotes - 5 - U.A.E. [24] 
S. doriae - 5 - Kuwait [14] 

S. doriae - 1 5 - Yemen [1 36] 
S. doriae - 6 - Israel [1 5] 
S. doriae -4 - Jordan [13] 
S. doriae - 14 - Yemen [136] ™ 
5. doriae - 3 - U.A.E. [12] 
S. doriae - 1 - Oman [11] 
S. doriae - 2 - Oman [6] 
S. slevini - 1 - Jordan [49] 

S.' slevini -i I Jordan [ 13 l 
S. slevini - 4 - Kuwait [50] 
S. steWn/ - 6 - U.A.E. [52] 
S. slevini - 5 - Qatar [51] 
r S. grandiceps - 3 - Jordan [17] 

!is. ff rand/ceps-2| Jordan [16] ^ 
l| S. grandiceps - 1 I V_y 
[nS.ffrafid/ceps- 9 1 Jordan [19] 
is. grand/eeps -10l L J 
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Figure 2 Bl tree of the genus Stenodactylus inferred using 12S, 16S mtDNA and c-mos, RAG-2 nuclear gene fragments. Black circles on 
the nodes indicate posterior probability values above 0.95 in the Bayesian Inference analysis. Numbers next to the nodes indicate bootstrap 
support of the Maximum Likelihood analysis (only values above 70 are shown). Ages of the nodes estimated with BEAST are indicated with an 
arrow, with the corresponding age range in brackets. The tree was rooted using Hemidactylus frenatus. Numbers in square brackets next to 
specimens code refer to localities in Figure 1. Information on the samples included is shown in Additional file 1: Table S1. Species' pictures were 
not submitted to precise relative scaling. 



monophyletic; (2) Tropiocolotes was forced monophyletic; 
and (3) Stenodactylus + Pseudoceramodactylus were forced 
monophyletic and Tropiocolotes was forced monophy- 
letic on the same constraint tree. The resulting con- 
strained topologies were compared to our optimal topology 



from Figure 2 under both ML and Bayesian frameworks 
(see Table 2). The results of the topological tests indicate 
that our dataset cannot reject the alternative hypothesis 
of monophyly of Stenodactylus + Pseudoceramodactylus 
(AU:0.461, SH:0.839, BF:0.647), monophyly of Tropiocolotes 
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(AU:0.161, SH:0.495, BF:-0.530) or both concurrently 
(AU:0.153, SH:0.492, BF:1.589). 

Within Stenodactylus, three well supported clades are 
revealed (see Figure 2): (i) clade A, formed by the Arabian 
species S. pulcher Anderson, 1896 [36], S. arabicus and 
the divergent lineage S. cf. arabicus, (ii) clade B, that 
includes five Arabian species {S. leptocosymbotes Leviton 
and Anderson, 1967 [83], S. doriae, S. slevini Haas, 1957 
[37], S. grandiceps Haas, 1952 [84] and S. affinis (Murray, 
1884) [85] grouped in 2 sub-clades, and (iii) clade C, 
formed by the four African species (S. petrii, S. stenurus 
Werner, 1899 [86], 5. mauritanicus Guichenot, 1850 [87] 
and S. sthenodactylus (Lichtenstein, 1823) and the south- 
west Arabian endemic 5. yemenensis Arnold, 1980 [31]. 

Clade A is sister to the remaining species of the genus 
and includes the two morphologically similar but highly 
divergent species S. pulcher and S. arabicus (p-distance 
12S: 12.5% and 16S: 14.5%) (Additional file 3: Table S2a). 
Genetic variability within S. arabicus is very high and 
includes two reciprocally monophyletic deep lineages 
(^-distance 12S: 7.7% and 16S: 5.0%) (Additional file 3: 
Table S2c), one of them restricted to the Sharqiya Sands 
(formerly Wahiba Sands) in Oman, hereafter referred to 
as S. cf. arabicus, and the other one covering the rest of 
the distribution range of the species. Network analysis of 
the nuclear gene fragments c-mos and RAG-2 shows that 
for the former all alleles are unique for each lineage and 
all but one for the latter (Figure 3). 

Clade B is well supported and groups S. doriae and 
S. leptocosymbotes in sub-clade Bl, while S. slevini, S. 
grandiceps and S. affinis in B2. Phylogenetic relation- 
ships are not completely resolved in the latter. Genetic 
distances between these five species are among the lowest 
in the genus (Additional file 3: Table S2a). Nuclear net- 
work analyses (Figure 3) reveal only unique alleles in 
the c-mos gene fragment for all five species, while there is 
some allele sharing in RAG-2 between S. doriae and 
S. leptocosymbotes. 



Finally, clade C consists of three sub-clades, two Afri- 
can and one Arabian. The North African sub-clade CI 
braches first, and the Arabian S. yemenensis is sister to 
sub-clade C3 formed by the two North African species 
5. mauritanicus and S. sthenodactylus, making the group 
of North African Stenodactylus species paraphyletic. 
Topological constraint analyses indicate that the alterna- 
tive hypothesis of monophyly of the North African 
species is rejected by the AU and BF tests (AU:0.029, 
SH:0.123, BF:7.221) (Table 2). 

In sub-clade CI, 5. stenurus is nested within S. petrii, 
rendering the latter paraphyletic. The results of the 
topological constraint analysis in which S. petrii was 
forced monophyletic show that this hypothesis is rejected 
by both AU and BF tests (AU:0.036, SH0.210, BF:2.578) 
(Table 2). Network analysis shows that S. stenurus lacks 
unique alleles in both nuclear markers (Figure 3). The 
level of intraspecific genetic variability within S. petrii 
(Additional file 3: Table S2b) is very high: the uncorrected 
/7-distances between specimens from Egypt and Israel, and 
the remaining S. petrii specimens sampled for this 
study is 7.2% and 6.0% for the 12S and 16S mitochon- 
drial markers, respectively (Additional file 3: Table S2c). 
Nuclear networks indicate that all six c-mos and four out 
of six RAG-2 alleles investigated are unique to this former 
lineage of S. petrii (Figure 3). 

In sub-clade C3, the two North African species S. sthe- 
nodactylus and S. mauritanicus are reciprocally mono- 
phyletic and highly divergent (/^-distance 12S: 10.9% and 
16S: 7.2%) (Additional file 3: Table S2a). The former is 
highly variable (p-distance: 12S 4.7% and 16S 3.2%) 
(Additional file 3: Table S2b) and presents three deep 
lineages that follow a clear geographical pattern (Figures 1 
and 2), grouping animals from: 1.- northern Egypt, Israel 
and Jordan; 2.- south, southeast Egypt and Kenya; 3.- all 
the animals from Libya, Algeria, Tunisia, Western Sahara 
and Mauritania, although a single specimen from NE 
Egypt (loc. 127 in Figure 1, Siwa Oasis) is also part of this 



Table 2 Statistical support for alternative hypotheses on Stenodactylus phylogeny 







ML framework 1 


Bayesian framework 2 


Tree 


-log likelihood 


AU P 


SH P 


HME 


log,o BF 


Unconstrained tree 


15180.095955 






-15175.4907 




Monophyly of Stenodactylus+Pseudoceramodactylus 


15180.877129 


0.461 


0.839 


-15175.0633 


0.647 


Monophyly of Tropiocolotes 


15183.765511 


0.161 


0.495 


-15175.6215 


-0.530 


Monophyly of Stenodactylus+Pseudoceramodactylus and Tropiocolotes 


15183.696084 


0.153 


0.492 


-15177.9132 


1.589 


Monophyly of African species 


15192.711115 


0.029 


0.123 


-15190.3457 


7.221 


Monophyly of S. petrii 


15189.220967 


0.036 


0.210 


-15181.4116 


2.578 



All topological tests are done versus the unconstrained (best) tree. Values in bold indicate statistically significant results. 

1 ML: Maximum likelihood; AU: Approximately Unbiased test (Shimodaira, 2002); SH: Shimodaira & Hasegawa (1999) test. P<0.05 suggests that the two solutions 
are significantly different. 

2 HME: The harmonic mean of sampled likelihoods as estimated by Tracer. BF: Bayes Factor. A log-| 0 Bayes factor > 2 indicates decisive evidence for statistically 
significant difference between solutions. 
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latter clade. The results of the network analyses also show 
a differentiation between these three lineages. In c-mos 
(Figure 3), in the first lineage 7 out of 10 alleles are 
unique, in the second lineage 2 out of 6 and in the third 
lineage 9 out of 16, while in RAG-2, 3 out of 10, 2 out of 8 
and 10 out of 12 are unique, respectively (Figure 3). 
Genetic variability within 5. mauritanicus is slightly 
higher than in S. sthenodactylus (^-distance: 12S 4.7% 
and 16S 4.3%) (Additional file 3: Table S2b) and six differ- 
ent mitochondrial lineages with geographic structure are 
found: 1.- easternmost part of Libya and Egypt; 2.- central 
Libya; 3.- Tunisia; 4.- Northern Morocco; 5.- two very 
divergent samples from southeastern Morocco; and 6.- all 
the southern Morocco plus Western Sahara samples (see 
Figure 1 and Additional file 2: Figure SI). 

Estimation of divergence times 

Convergence was confirmed examining the likelihood 
and posterior trace plots of the two runs with Tracer 
v.1.5. Effective sample sizes of the parameters were 
above 200, indicating a good representation of independ- 
ent samples in the posterior. The estimated divergence 
times are illustrated in Figure 2 and the chronogram can 
be seen in Supplementary Figure 2 (Additional file 4: 
Figure S2). Diversification within Stenodactylus initiated 
29.5 Ma ago (95% HPD: 20.7-39.2). In clade A, the split 
between S. pulcher and S. arabicus is dated back to 
17.3 Ma (95% HPD: 11.3-23.6). The separation between 
the ancestors of clades B and C dates back to 21.8 Ma 
(95% HPD: 15.4-29.1), while diversification within these 
two clades started 11.0 Ma (95% HPD: 7.4-15.0) and 
19.3 Ma (95% HPD: 13.1-25.5) ago, respectively. 

Posterior mean rates for the 12S and 16S mitochon- 
drial gene fragments were estimated at 0.00701 and 
0.00642 substitutions per lineage per million years, re- 
spectively (or divergence rate: 1.402% and 1.284%). The 
posterior rates for the nuclear fragments, c-mos and 
RAG-2, were 0.00052 and 0.00060 respectively, more 
than 10 times lower than the mitochondrial ones. The 
12S mitochondrial rate concords extremely well with the 
average rate for the same mitochondrial gene for three 
Canary Island reptile groups (Gallotia, Tarentola and 
Chalcides; 0.00755 for the 12S gene) as estimated by 
Carranza and Arnold (2012) [44]. 

Ancestral area reconstruction 

Reconstruction of the ancestral areas of Stenodactylus 
species was done in a parsimony framework based on 
the topology of the phylogeny presented in Figure 2. The 
analysis indicates that the reconstruction of the area for 
some of the ancestors is equivocal (see Figure 4). These 
are the common ancestor of clade C, formed by all North 
African species and S. yemenensis, and the ancestor 
of the latter and the sister species S. sthenodactylus I 



S. mauritanicus. Reconstructions using accelerated 
transformation (ACCTRAN) or delayed transform- 
ation (DELTRAN) optimizations support an identical 
number of events involving Arabia and Africa, but 
the direction of events is different. ML-based reconstruc- 
tion, considering branch-length information, with the 
best- fit Markov k-state 1 -parameter model also provided 
results with fairly similar probabilities for the two states in 
the aforementioned nodes (Figure 4). 

Discussion 

This constitutes the first phylogenetic study using a 
complete sampling of Stenodactylus taxa and including 
207 specimens from across the entire distribution range 
of North Africa and Arabia (Figure 1). This has enabled 
a robust phylogenetic reconstruction (see Figure 2 and 
Additional file 2: Figure SI), the uncovering of intraspeci- 
fic diversity and, in some cases, the unveiling of interesting 
distribution patterns (see below). The phylogenetic results 
show a high level of support in most of the nodes and a 
striking agreement with the phylogenetic analyses of 
Stenodactylus by Arnold (1980) [31], based on morpho- 
logical data, increasing our confidence that the recovered 
topology represents the true evolutionary history of 
the genus. 

Monophyly of Stenodactylus 

Despite the general concordance between morpho- 
logical and phylogenetic conclusions, one important 
discrepancy is observed: while morphology supports 
the inclusion of P. khobarensis in the genus Stenodactylus, 
the results of our molecular analyses indicate that 
Pseudoceramodactylus and Stenodactylus are not even 
sister genera (Figure 2). Kluge (1967) [39] transferred 
P. khobarensis to the genus Stenodactylus based on a 
"large number of external (meristic and mensural) and in- 
ternal morphological similarities", including relevant cha- 
racters like the phalangeal reduction to a formula of 
2.3.3.4.3 on both fore and hind limbs and a very high scleral 
ossicle number (20-28). Arnold (1980) [31], despite poin- 
ting out some unique scale characters of P. khobarensis, 
retained it in Stenodactylus and considered the scalation 
characters as "convincing pointers to holophyly". However, 
according to a recent molecular analysis of the group by 
Fujita and Papenfuss (2011) [40] based on independent 
samples and sequences of different mitochondrial and nu- 
clear regions, two representatives of Tropiocolotes branched 
between P. khobarensis and the six species of Stenodactylus 
included in the analysis (see Figure 1 of [40]). In order to 
deal with the non-monophyly of Stenodactylus, the genus 
Pseudoceramodactylus was resurrected. This pattern is 
repeated and further investigated in our study, with a 
complete taxon sampling of Stenodactylus and the inclu- 
sion of a greater number of representatives of Tropiocolotes, 
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Figure 4 Ancestral area reconstruction. The tree figure illustrates the parsimony reconstruction, while numbers above and below nodes 
correspond to ML probabilities for character states. Black and white colors correspond to Africa and Arabia respectively, and grey color indicates 
equivocal nodes. 



resulting in the splitting of the latter genus into two 
groups, a surprising but not strongly supported, albeit 
consistent, result. 

We performed a series of constraint analyses in which 
Stenodactylus and Pseudoceramodactylus were forced 
to form a monophyletic group. Results clearly show 
that our dataset cannot reject the alternative hypothesis of 
a monophyletic Stenodactylus + Pseudoceramodactylus 
group (Table 2). In order to further investigate this, 
the dataset of Fujita and Papenfuss (2011) [40] was 
subjected to the same ML topological tests, but also 
could not reject the alternative hypothesis of mono- 
phyly of Stenodactylus + Pseudoceramodactylus (AU 
P = 0.074; SH P = 0.092). In view of the confusing molecu- 
lar evidence and taking into account the morphological 
data, we think that the resurrection of Pseudoceramodac- 
tylus was precipitated, but in the meanwhile, this change 
accommodates for both the paraphyly reported by Fujita 
and Papenfuss and confirmed here, and the hypothesis of 



monophyly of Stenodactylus + Pseudoceramodactylus. We 
recommend not performing any further changes at the 
generic level before an in-depth revision clarifies the evo- 
lutionary relationships between the genera Stenodactylus, 
Pseudoceramodactylus and Tropiocolotes. 

Systematics and evolution 

The well-supported clade A is formed by the morpho- 
logically similar S. pulcher, S. arabicus and the lineage 
S. cf. arabicus and, according to the inferred dates, 
the split between the former and the two latter spe- 
cies dates back to approximately 17 Ma ago (95% 
HPD: 11.3-23.6) (Figure 2). On the one hand, variabil- 
ity within S. pulcher is very low, probably as a result 
of the two specimens analyzed being from very close 
localities. On the other hand, the S. cf. arabicus 
lineage from the Sharqiya Sands (formerly Wahiba 
Sands), as already highlighted by Fujita and Papenfuss 
(2011) [40], is genetically very distinct from all other 
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populations of S. arabicus included in our study in 
both mitochondrial and nuclear markers (Additional 
file 3: Table S2c and Figure 2), where almost all alleles 
are lineage-specific (see Results and Figure 3). This sup- 
ports the idea that the Sharqiya Sands are isolated and 
surrounded by some areas of unsuitable habitat for sand 
dune specialists like this species [88-90]. Further morpho- 
logical and molecular studies including more specimens 
from putative contact zones and faster nuclear markers 
are expected to give S. cf. arabicus formal recognition. 

Clade B is well-supported (ML 100%, BI 1.0) and was 
also recovered by the morphological analysis of Arnold 
(1980) [34]. Stenodactylus doriae and S. leptocosymbotes 
are reciprocally monophyletic and form the relatively 
well-supported sub-clade Bl (Figure 2). Our molecular 
results agree with the results of the morphological ana- 
lysis by Arnold (1980) [31], who also recovered the two 
species as sister taxa based on three synapomorphies. 
The two species diverged approximately 7.0 Ma ago (95% 
HPD: 4.2-10.1) (Figure 2) and, like the two North African 
sister species S. sthenodactylus and S. mauritanicus, they 
are ecologically distinct. Stenodactylus leptocosymbotes is 
an arid-adapted species that lives on relatively hard, 
although usually sandy, substrates being replaced by its 
sister species, 5. doriae, on soft, wind-blown sand [34,91]. 
Thanks to its morphological and physiological adapta- 
tions, the latter is able to live in hyper-arid sand dune 
environments like for example the Eastern Rub al Khali 
[92], one of the largest and driest sand deserts in the 
world [93]. Given the clear morphological and ecological 
differences between these two species and the apparent 
absence of morphologically intermediate individuals 
[31,34], it seems reasonable to deduce that allele sharing 
in RAG-2 (see Results), which is limited to the ancestral 
allele, is the result of incomplete lineage sorting rather 
than ongoing gene flow between the two species. Varia- 
bility within S. leptocosymbotes is rather low (Additional 
file 3: Table S2b) and the number of samples included per- 
mit to observe only moderate geographical structuring 
(Figures 1 and 2, Additional file 2: Figure SI). In contrast, 
S. doriae, shows a higher level of genetic differentiation, 
with the Sharqiya Sands lineage being quite divergent 
(Additional file 3: Table S2c and Figure 2), as already men- 
tioned by Fujita and Papenfuss (2011) [40]. 

Sister to sub-clade Bl is a group composed by S. slevini, 
S. grandiceps and S. affinis, for which support is relatively 
low (ML 62, BI = 0.95). The topology within this sub-clade 
differs from the morphological hypothesis of Arnold 
(1980) [34], which supported the following relationship: 
(S. grandiceps (S. affinis {S. slevini (S. leptocosymbotes, 
S. doriae)))). Stenodactylus slevini is the only member 
of the group with two divergent lineages, one limited 
to Jordan and the other with representatives from 
East Arabia. Although the divergence based on mi- 



tochondrial data is clear (Additional file 2: Figure SI), 
there is no supporting nuclear data available (Figure 3), 
and no obvious morphological differences (pers. obs.). 
With the only exception of the soft wind-blown sand spe- 
cialist S. doriae, all remaining representatives of clade B 
plus two other species, the African S. sthenodactylus and 
the Arabian S. yemenensis, appear to occupy rather similar 
spatial niches. These six species are adapted to living on 
relatively hard ground, coarse sandy planes, large wadis 
and sandy substrates and, based on their head dimensions, 
probably feed on similar-sized prey [31,32,34]. As a conse- 
quence of that, these species rarely coexist and have 
largely allopatric distribution ranges, while in places where 
they coincide they are not syntopic [31,33,34]. The ana- 
lysis of the nuclear allele networks (Figure 3) indicate that 
the morphologically and ecologically similar and phylo- 
genetically closely related S. leptocosymbotes, S. slevini, 
S. grandiceps and S. affinis do not share a single allele 
in the c-mos and RAG-2 genes analyzed, even though the 
results of the calibration analyses suggest that S. grandiceps 
and S. affinis diverged later (6.7 Ma ago; 95% HPD: 4.1- 
9.3) than other lineages for which extensive allele 
sharing in the RAG-2 has been identified (S. doriae 
and S. leptocosymbotes; see above and Results). These 
differences of the level of lineage sorting in some of 
the morphologically well-recognized species may also 
be the result of differences in effective population 
sizes, which affect the lineage coalescence time [94]. 

In sub-clade CI, S. petrii is grouped together with the 
North African endemic S. stenurus that branches inside 
it (Figure 2). As a result, S. petrii is paraphyletic and 
constitutes the only exception among the otherwise 
monophyletic Stenodactylus species. The results of 
the topological tests (Table 2) indicate that our data- 
set most probably rejects the monophyly of this species 
(AU:0.036, SH:0.210, BF:2.578). Stenodactylus stenurus 
was described by Werner (1899) [86] and synonymized 
ten years later by the same author [95]. It remained in syn- 
onymy until Kratochvil et al. (2001) [96] recognized it as a 
valid species, based on a multivariate analysis of several 
metric and scalation characters. It is noteworthy that the 
representative of S. stenurus included in our analysis is 
one of the specimens used by Kratochvil et al. (2001) [96] 
in their study. 

The highly divergent lineage that includes specimens 
from Egypt and Israel (see Results) is estimated to have 
split from specimens further west in Algeria, Morocco, 
Western Sahara and Mauritania approximately 6.1 Ma 
ago (95% HPD: 3.9-8.6) (Figure 2). In fact, the northern 
Sinai populations of S. petrii have been reported to be 
morphologically distinct and, as a result of that, were 
considered a different species (S. elimensis) by Barbour 
(1914) [97], now under the synonymy of S. petrii [31,98]. 
Yet, specimens from this area included in our analyses 
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do not present considerable genetic differences with the 
rest of the Egyptian and Israeli specimens (Figures 1 and 
2, Additional file 2: Figure SI). It should be pointed out 
that the type locality of S. petrii is Egypt and, thus, this 
lineage represents the 'true' S. petrii. The pattern in 
the nuclear genes, with numerous unique alleles for 
this lineage (Figure 3), contrasts with the situation in 
S. stenurus that lacks unique alleles. This suggests 
that further analyses and a thorough taxonomic revi- 
sion including more samples of S. petrii, especially 
from not sampled areas of Algeria and Libya, and 
mainly S. stenurus will be necessary in order to evalu- 
ate the status of the populations assigned to the two 
species. With this evidence it will be possible to dif- 
ferentiate between a single species with high genetic 
variability (petrii), two species (petrii in the East and 
stenurus in the West) or three species, if stenurus 
proves to be distinct from the more western forms. 

The two North African species of sub-clade C3, 
S. sthenodactylus and S. mauritanicus, are shown to be 
reciprocally monophyletic and highly divergent (Additional 
file 3: Table S2a), while their separation dates back to ap- 
proximately 10.0 Ma (95% HPD: 6.6-13.7) (Figure 2). These 
results help to clarify the status of these two taxonomically 
controversial taxa that were treated as two different subspe- 
cies by Loveridge (1947) [99] and Sindaco and Jeremcenko 
(2008) [33], as the same monotypic species by Arnold 
(1980) [31] and that were finally considered as full species 
by Baha el Din (2006) [98], who found them in sympatry at 
particular localities in northern Egypt. As observed by Baha 
el Din (2006) [98], although these two sister species can be 
morphologically similar and share similar habits, they are 
ecologically different. Stenodactylus mauritanicus is 
restricted to fairly mesic coastal semi-desert under the in- 
fluence of the Mediterranean (see Figure 1), where it inha- 
bits flat rock-strewn sand and gravel plains with fairly good 
vegetation cover. On the contrary, S. sthenodactylus inha- 
bits areas of the Sahara that are far more arid and inhospit- 
able than the ones of its sister species (see Figure 1), being 
the only vertebrate to be readily found in some parts of the 
Western Desert of Egypt [98]. It prefers gravelly and coarse 
sandy plains and large wadis and, although the species is 
typical of hard coarse substrates, it sometimes penetrates 
some dune areas [98]. 

The distributions of these two species, as introduced by 
the present study, give insights into the controversial taxo- 
nomic status and frequent misidentification of the two 
forms [99]. Our analysis concludes that S. sthenodactylus 
extends west from the Middle East and Egypt, previously 
thought to be its eastern limit, across the Sahara and into 
Mauritania (Figure 1). Stenodactylus mauritanicus is con- 
firmed to be present in Egypt [98] and has a wide, almost 
continuous distribution roughly along the northern margin 
of the Sahara desert. The two species are found in 



sympatry or in close proximity in Egypt and coastal 
Mauritania, yet retain distinct mtDNA lineages and exhibit 
only limited allele sharing in the nuclear markers, most of 
which is due to sharing of ancestral alleles and hence is 
likely to represent incomplete lineage sorting (see 
Figure 3). 

Stenodactylus sthenodactylus presents high variability, 
both at genetic (see Results) and morphological levels 
[31]. Its three deep lineages are estimated to have 
diverged approximately 4.8 Ma ago (95% HPD: 2.8-6.9) 
(Figure 2). According to Baha el Din (2006) [98], some 
morphological characters appear to correlate with envi- 
ronmental factors, with populations from hyper-arid places 
showing a very slender body, less contrasting pattern and 
tubular nostrils, while populations from more mesic areas 
being usually more robust, with thick limbs, big heads and 
marked pattern [31,36,98]. The populations from coastal 
regions in southeast Egypt are especially distinct and, 
according to Baha el Din (2006) [98], they resemble speci- 
mens of 5. s. zavattarii from Kenya, which Loveridge 
(1957) [100] synonymized with 5. sthenodactylus. Two spe- 
cimens of this form were included in our phylogenetic ana- 
lyses (see Figure 2 and Additional file 2: Figure SI), 
and indeed they belong to a clade with samples from 
south and southeast Egypt. These results suggest that 
some of the morphological variability between popula- 
tions of S. sthenodactylus may also be supported by mo- 
lecular data. A nomenclatural revision of North African 
Stenodactylus (work in progress) is essential for stability 
before any changes are performed, while further work 
focused on the contact zones between the three lineages 
and combining detailed morphological analyses with add- 
itional nuclear data is needed in order to determine if they 
deserve formal recognition. 

On the other hand, the high genetic variability within S. 
mauritanicus (Figure 2 and Additional file 3: Table S2b) 
does not seem to correlate with differences in morphology. 
This species is fairly uniform morphologically, with 
populations from the West being a bit larger than 
Egyptian ones but generally maintaining the same propor- 
tions, pattern and scalation across most of its distribution 
range [98]. Nevertheless, the intra-specific divergence is 
estimated to date back to 6.6 Ma ago (95% HPD: 4.0-9.5) 
and the six mitochondrial lineages present a clear geo- 
graphical pattern (Figure 1 and Additional file 2: Figure SI). 
The relationship between these lineages, however, is not 
clear and neither is any structure observed in the nuclear 
alleles (Figure 3), both facts being mirrored in the low- 
supported nodes of the concatenated phylogeny (Figure 2). 

Origin, biogeography and diversification of Stenodactylus 

Reconstruction of ancestral areas with both parsimony 
and ML methods (Figure 4) suggests that the genus 
Stenodactylus originated in Arabia approximately 
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30 Ma ago (95% HPD: 20.7-39.2) (Figure 2), a time of 
high geological instability as a result of the onset of 
major seismic and volcanic events in the general area 
of Ethiopia, northeast Sudan and southwest Yemen 
[101]. These major volcanic and tectonic events, cen- 
tered over the Afar region, marked the onset of the 
formation of some of the most relevant and complex 
physiographical features in the contact zone between 
Africa and Arabia, like the Gulf of Aden, the Red Sea 
and the elevation of the Afro-Arabian rift-flanks to 
heights above 3600 m [1,101,102]. 

The tempo and mode of the deep splits in Stenodactylus 
bear a striking resemblance to the basal splits that oc- 
curred in the African-Eurasian snake genus Echis [13,103], 
which suggests a common biogeographical pattern 
for both groups. The distribution of the members of 
Arabian clade B (S. doriae, S. leptocosymbotes, S. slevini, 
S. grandiceps, S. affinis) and the mainly African clade C 
(S. petrii, S. stenurus, S. yemenensis, S. mauritanicus, 
S. sthenodactylus) (Figures 1 and 2) extend primarily 
on the opposite sides of the Red Sea, mimicking the situ- 
ation of the sister taxa E. coloratus (mainly Arabian) and 
E. pyramidum (mainly African). The split between these 
two Stenodactylus groups dates back to 21.8 Ma ago (95% 
HPD: 15.4-29.1) (Figure 2), which roughly coincides with 
the split between E. coloratus and E. pyramidum calcu- 
lated at approximately 19.4 Ma ago. The dates of these 
phylogenetic events follow a well-studied phase of volcan- 
ism and strong rifting initiated at approximately 24 Ma 
ago, that appeared in an almost synchronous way through- 
out the entire Red Sea [1]. Therefore, it is possible that the 
formation of the Red Sea acted as a vicariant event separa- 
ting the aforementioned clades of Stenodactylus, as also 
suggested by Pook et al. (2009) [13] for the genus Echis. 
The agamid lizards of the genus Uromastyx [25] is yet an- 
other group that could have been affected by such an 
event, although in this case the split between the 
Arabian and African clades seems to have happened 
later, at 11-15 Ma ago. Amer and Kumazawa (2005) [25] 
attributed this split to a dispersal event from Arabia into 
North Africa, coinciding with climatic changes towards 
aridity in this latter area, rather than to vicariance. How- 
ever, since earlier dates had also been calculated for the 
split between African and Arabian Uromastyx that coin- 
cide with the inferred dates for Stenodactylus and Echis 
(18 Ma ago; [104]), a reassessment of the calibration dates 
of Uromastyx using relaxed clock methods like the ones 
applied by Pook et al. (2009) [13] and in the present study 
seems necessary (work in progress). 

The split between the Arabian S. yemenensis and the an- 
cestor of the African S. mauritanicus and S. sthenodactylus 
on either sides of the Red Sea also parallels the 
splits between Arabian and African sister clades of the 
E. pyramidum complex [13] and Uromastyx ocellata and 



U. ornata [25]. Although the divergence time estimate for 
the Stenodactylus members (15.4 Ma ago (95% HPD: 
10.5-20.8), Figure 2) predates the ones of the other two 
groups by almost 7 Ma, the split between African and 
Arabian lineages might be explained by the complex geo- 
logy of the Red Sea. Several recurrent episodes during the 
Miocene caused the desiccation and refilling of this tec- 
tonically active rifting area [1,105] and provoked the seve- 
ring of the land bridges that had existed after the initial 
formation of the Red Sea in the early Miocene. So, the 
separation between S. yemenensis and the ancestor of 
S. mauritanicus and S. sthenodactylus was probably 
also the result of vicariance, similarly to Echis and 
Uromastyx. After this event, 5. yemenensis would have 
remained isolated at the coastal side of the southern 
Arabian highlands (Figures 1 and 2). 

In Arabia, an example of a similar biogeographical pat- 
tern caused by a different biogeographical process is the 
case of the ecologically similar sister species of clade A, 
S. pulcher and 5. arabicus (including S. cf. arabicus), 
which, according to the results (Figure 2) and the geo- 
logical data available, are hypothesized to result from 
vicariance caused by the uplift of the Yemen Mountains 
approximately 18 Ma ago [1,101,102]. The splits within 
clade B, however, seem more difficult to interpret, as little 
information is available on the geological and climatic his- 
tory of the interior of Arabia. A general pattern could be 
proposed with a first North-South split between the 
ancestors of S. doriae, S. leptocosymbotes and S. slevini, 
S. grandiceps, S. affinis, respectively, followed by the pos- 
terior range expansion of some of these species. Interest- 
ingly, in Arabia, even though evidence exists for an increase 
in aridification [106], it has been hypothesized that at the 
same time an important river system, as evidenced by the 
fluvial sediments, could characterize the interior of the pe- 
ninsula [93,107]. Such dynamic scenery could be respon- 
sible for the rapid diversification within clade B, having 
caused fragmentation of the distribution range of the ances- 
tor^) and the different lineages to split allopatrically. 

The onset of diversification in clade B coincides in 
time with the split between the African S. mauritanicus 
and S. sthenodactylus in sub-clade C3 (Figure 2). These 
speciation events match very closely the estimates of the 
formation, in the late Miocene, of a major east- Antarctic 
ice sheet with its associated polar cooling, which triggered 
the aridification of mid-latitude continental regions and a 
shift in North Africa from forest to dry open woodlands 
and savannahs [4,20,108]. The two North African forms, 
S. mauritanicus and S. sthenodactylus, seem to have 
diverged in ecological niche, with one form adapted to 
mesic environments and the other occupying much dryer 
areas, respectively. It has been proposed that the gradual 
increase in aridity that took place in northern Africa 
during the late Miocene accelerated the diversification 
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process in reptile faunas [21]. The estimated diver- 
gence times of the North African Stenodactylus seem 
to corroborate a common emerging pattern among 
European biota, according to which the speciation 
events in many reptile and amphibian groups do not 
coincide with the accentuated environmental instability 
during the Pleistocene, but rather date back into the 
Miocene and proceeding through the Quaternary, when 
many species and populations originated [109,110]. 

It has been suggested that 18 Ma ago, Africa con- 
nected with Eurasia through the closure of the Eastern 
Mediterranean seaway (the Gomphotherium land bridge) 
[15]. This land bridge later became disconnected temporar- 
ily but it has been continuously present since approximately 
15 Ma ago. It is interesting to notice that, despite the exist- 
ence of a continuous passage between Arabia and Eurasia, 
our phylogeny suggests that colonization of Eurasia by 
members of the genus Stenodactylus occurred much later 
and was very restricted geographically. In fact, only two 
Stenodactylus species extend their ranges into Eurasia 
(S. affinis and S. doriae). From these two, only samples 
of S. affinis from Eurasia (Iran) were available, while for 
the other species a specimen from neighboring Kuwait was 
included. In both species, however, the low intraspecific 
genetic variability suggests that the colonization of Eurasia 
was a very recent event (Figure 2 and Additional file 3: 
Table S2b). One possible explanation of this biogeographical 
pattern may be the existence of ecologically and morpho- 
logically very similar forms in Iran like Crossobamon 
(formerly a member of Stenodactylus [39]) and Agamura, 
which may compete with Stenodactylus and therefore may 
have not allowed it to expand further outside the narrow 
coastal strip in southwestern Iran (Arnold, 1980). This situ- 
ation is completely different than the one in North Africa, 
where no ecological analogs to Stenodactylus seem to exist 
and therefore several of its species are found across an area 
of more than 10 million Km 2 [31,33,98,111,112]. 

Conclusions 

The analyses presented in this study, based on a multi- 
locus dataset that derives from a complete sampling of 
the 12 species of the genus Stenodactylus, reveal the 
existence of three clades with deep divergences within 
Stenodactylus and high intraspecific variability in some 
species, while the estimation of divergence times allows for 
biogeographical interpretations. The geckos Stenodactylus 
originated in Arabia 30 Ma ago. In clade A, the split be- 
tween the two species is hypothesized to have resulted 
from vicariance caused by the uplift of the Yemen Moun- 
tains approximately 18 Ma ago. Stenodactylus cf. ambicus 
from the Sharqiya Sands constitutes a genetically and mor- 
phologically distinct lineage. In clade B, rapid diversifica- 
tion seems to relate to climatic and geological instability in 
the late Miocene, but this hinders the reconstruction of 



robust phylogenetic relationships between some species. 
The Sharqiya Sands host yet another divergent lineage, that 
of the species S. doriae. In clade C, the split between 
S. yemenensis and sub-clade C3 is hypothesized to re- 
late to the recurrent episodes of the desiccation and 
refilling of the Red Sea, during the Miocene. An interest- 
ing distribution pattern is revealed for the sister species 
S. sthenodactylus and S. mauritanicus, differing greatly 
from what was previously thought. Several speciation 
events in Stenodactylus are estimated to date back to 
the late Miocene, indicating that this was an impor- 
tant period for reptile diversification in this area. The 
split between clades B and C is attributed to the 
opening of the Red Sea in the Upper Miocene, acting 
as a vicariant agent. On the other hand, the forma- 
tion of the connection between Africa and Eurasia 
seems to have had little effect on Stenodactylus, pro- 
bably because of the existence of ecological analogs. 
On a taxonomic level, further studies are expected to re- 
solve the systematics of the S. petrii - S. stenurus complex. 
Validity of the specific status of S. mauritanicus is con- 
firmed with mitochondrial and nuclear data. Overall, this 
work unveils the evolutionary history of Stenodactylus 
geckos and highlights their use as a model in the study of 
the faunal interchanges between North Africa and Arabia 
and the evolutionary processes in these arid areas. 
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